Conversion of free magnetic polaron into vortex lattice in diluted magnetic 
semiconductors in quantizing magnetic fields. 
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We show that in strong (quantizing) magnetic fields "ordinary" free magnetic polaron in diluted 
magnetic semiconductors (of type A 2 MeB 6 , where Me = Fe, Mn) exists in the form of vortex lattice 
quite similar to that in type II superconductors (Abricosov vortex lattice). The region of external 
parameters (like external magnetic field and temperature), where such lattice exists, is determined 
from the condition that lattice dimension is less or equal to polaron localization radius. 
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Diluted magnetic semiconductors (DMS) can be re- 
garded as an ensemble of magnetic impurities Fe, Mn, 
embedded in the lattice of host A 2 B 6 semiconductor. The 
concentration of impurities is such that direct exchange 
interaction between impurities does not occur so that 
aforemetioned substances is in the paramagnetic phase. 

Last decade considerable experimantal and theoretical 
efforts have been spared for investigation of self-localized 
states of carriers in DMS (see, e.g. Such self- 

localized state of conduction electron (hole) is called free 
magnetic polaron (MP) and can be regarded as a region 
(of polaron radius p) of correlated spins of magnetic ions, 
where the interaction between spins of magnetic impuri- 
ties is mediated by the electron (hole). 

However, existing theories considered the properties of 

such MPs in magnetic fields (such that pu = *J jff < p, 
p is MP localization radius, H is external magnetic field), 
where diamagnetic effects (Landau quantization) do not 
manifest themselves (see (H)). It turns out, however, that 
shape and physical properties of MP may change drasti- 
cally if we include aforementioned diamagnetic contribu- 
tion into consideration. 

In this paper we show that MP in quantizing magnetic 
fields exists as a lattice of vortices, very similar to Abri- 
cosov vortex lattice in type II superconductors ||. We 
have shown that formation of the vortex lattice is ener- 
getically favourable as compared to the single MP. Our 
analysis show that general scenario of vortex lattice for- 
mation is as follows. When magnetic field exceeds some 
threshold value 

H cl = — 2 , (1) 
ep A 

the initial spherically symmetric "paramagnetic" MP 



(see H for detailed investigation of this MP properties; p 
is a function of magnetic ions concentration, temperature 
and external magnetic field) looses its sprerical symme- 
try and gains cylindrical one. This "diamagnetic" MP is 
localized in a plane perpendicular to magnetic field di- 
rection (i.e. in the plane of Landau quantization). At 
further increase of magnetic field this MP splits into sev- 
eral vortices, that tend to organize into vortex lattice 
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We can draw an "analogy" between such "diamag- 
netic" MP and vortex in type II superconductor. While 
magnetic field penetrates only to the vortex core in super- 
conductor, giving rise to normal (nonsuperconducting) 
phase there, the core region of MP consists of coherently 
oriented correlated spins of magnetic impurities, giving 
rise to (local) magnetically ordered phase. 

At very high magnetic fields H ~ H C 2 the "ferromag- 
netic" cores of the vortices merge so that magnetic impu- 
rities do not "need" the electron (hole) to interact with 
each other. In this region of magnetic fields the ampli- 
tude of MP wave function become small so that H C 2 can 
be determined by the usual procedure of linearization 0] 
of corresponding equation. 

In the present paper we give the estimations of H c i 
and H C 2 rather then detailed (standard) calculation of 
vortex lattice properties. 

Consider first the model of free MP. We start from the 
Hamiltonian of electron (hole) , coupled to subsystem of 
magnetic ions in DMS. These magnetic ions randomly 
substitute Ni cations at the sites Rj (j = 1,2, ...Ni) of 
a cubic crystal lattice. Under the assumption that cor- 
responding energy band is isotropic and nondegenerate, 
we write the Hamiltonian of this system in the effective 
mass approximation 



H = 



2m* 



Ni 

g e p B HS e + ]T [-JQS (f- 4) SiS e 



(2) 
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where p, m* , g e , S e and r are momentum operator, effec- 
tive mass, band g-factor, spin operator, and coordinate 
of electron (hole); gi and Sf are g-factor and spin op- 
erator of the magnetic ion localized at lattice site Rj, 
J is carrier-ion exchange interaction constant, Q is unit 
cell volume and is Bohr magneton. The first two 
terms in (|J) determine the kinetic and Zeeman energy of 
an electron. The next two terms determine the carrier- 



ion exchange interaction and the sum of magnetic ions 
Hamiltonians in the field H (see || for details). 

Choosing electron wave function as a product of coor- 
dinate and spin parts and averaging over random impu- 
rities positions in a mean field approximation (see Q for 
details), we obtain following energy functional for a free 
magnetic polaron. 



E„ = 




Pz--A z ) il> + - 

c / 2mj_ 



+ m J [Ei (gip B H + Qa) - E x (gi^ B H)] d 3 r, 



ip* Pi A±) ijj >d r + g e n B H(T 



(3) 



where a = ±i are the eigenvalues of the operator pro- 
jecting the electron spin onto H and 



'(f) = -Jft|V>(f)| 2 



(4) 



is the effective exchange field, E\ depends on the specific 
type of SMS magnetic ions paramagnetism (see below), 
p z = -ihJ^,p ± = -ih(e x J^+e y -^), A± = e x A x +e y A y . 

For SMS of type A 2 _ x Fe x B 6 with Van Vleck paramag- 
netism Ei has the form j8| 



Eiw(x) =Ae- [(Ae) 2 +4x 2 ] 



1/2 



(5) 



where Ae is the splitting of the spin-orbit multiplet of 
the magnetic ion (see || for details). For SMS of type 
A 2 _ x Mn x B 6 with orientational paramagnetism we have 
for Ei 



ElOR (X) - - lOg -r—^ 

P sinh (f) 



(G) 



where S is magnetic ion spin, (3 — {k B T) 

It is seen that (||) looks like Ginzburg-Landau (GL) 
functional for superconductors. The difference is more 
complicated nonlinearity. It is seen that expansion of 
Eqs (|J) and (0) in powers of \ip (r)\ 2 up to (f)^ put 
formal one-to-one correcpondence with GL functional. 

Variation of (||) with respect to ip* gives following 
Schrodinger equation for polaron wave function 



ip + m 



dEi 



i) = o. 



(7) 



Here we omit the unimportant dependence on z, since 
the electron motion is not quantized in this direction. 
Boundary conditions for MP are similar to those for GL 
superconducting order parameter and require that ip = 
at infinity. 



It is convenient to show the existense of vortex lattice 
solutions of corresponding differential equations from the 
side of H ~ H C 2, where \ip (f)| 2 is small. For that we lin- 
earize (|?]) (with respect to (||) and (|^)) in ip. This gives 



1 { 6 \ 

t; l-ifi— + -Hx) ip + Xip = 0, 

2m± \ tiy c J 

where we choose following gauge of vector potential 



A„ 



Hx. H = H 7 



(8) 



(9) 



wv 



Agifj, B H 



v /(Ae) 2 +4( 3 ^ BJ ff)2' 
X OB = -jnn l aSB s (Sf3g l fi B H), (10) 



Bs(x) is a Brillouin function for spin S: 
SB s (Sx) = f S + i J coth ( S + i J x 



1 , x 
— coth — . 
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It is well-known (see, e.g. jjj) that localized solutions of 
(||) exist if A equals to one of eigenvalues 

A=- (n+- )h . 

2 m_i_c 

Lowest eigenvalue determines H C 2 so that we obtain 
following equations for H C 2 for the case of Van Vleck 



Jxio- 



4gi/j, B H 



1 eH 

-h- 



V(A£) 2 +4(«?^ B tf) 2 2 m ± c 
and orientational paramagnetism 

<7JxiSB s {SPgi^ B H c2 oR) - ^ bHc2 ° r = o, 

m± 

/i = — , xi = mil, 

where mo is a free electron mass. 



(11a) 



(lib) 
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To have feeling of the numerical values of H C 2 for 
some specific DMS, we use two examples. First one 
is Zni_ x Fe x Se with Van Vleck paramagnetism of Fe. 
Its parameters are ]To| : m* = 0.14too; m h = 1-2too; 
J e = 0.22eV; J h = -1.6eV; Ae = l.SmeV; unit cell vol- 
ume O = 45. 6 A 3 . It is sufficient for our estimations to 
put m* h = mxe,h- 

Second one is Cdi_ x Mn x Te with orientational param- 
agnetism of Mn. Parameters of this solid solution are fol- 
lowing Q: m* = 0.096m ; to* = 0.48mo; J e = 0.22eV; 
Jh = — 0.88eV; lattice constant a = 6A8A; unit cell vol- 
ume Q = 68. 06 A 3 . 



It is seen from (11a) that this equation has solution 
only if J a < 0, i.e. for vortex lattice to occur for electron 
MP (J e > 0), we should have a = -1/2; for hole MP 
(7 = 1/2. We have from (11a) 



1 



1 



Ae 



4h 2 c2VV 



4 [igixi\Ja\ 



Ae ' 



(12) 



Since Ae/J « 1, the solution of (|T^) is at h >> 1. In 
this case we have from (n2|) 



H C 2vv = (1-7 • 10 4 Tesla) fjixi\J \, J = J/leV (13) 

It was shown earlier Q that in Zn!_ x Fe x Se electron can- 
not autolocalize. So, we give estimate of H c2 for hole 
only. From (|l^) (at typical value xi — 4%) we obtain, 
that for hole H^ 2VV » 1300 Tesla. 

The solution of Eq. ( |llh| ) exists at J a > 0. In dimen- 
sionless variables 



h, 



3k B T 



JS{S + l)/JLXl 



equation (lib) can be rewritten as 

3h c 20R 



lc20R 



= B, 



t{S+1) 

Its asymptotic values are as follows 




0.4 0.6 0.8 

3k B T/JS(S+1)nx 

FIG. 1. "Universal" temperature dependence of H C 2 for 
orientational DMS. 

The situation with H c \ is quite similar to that in type 
II superconductors. In this case MP plays a role of single 
vortex. Our analysis show that at H < H c \ MP exists in 
3D spherically symmetrical form while at H > H c i it has 
2D character being localized in the plane perpendicular 
to the direction of the external magnetic field. So, for H c i 
estimation it is sufficient to consider spherically symmet- 
rical MP, substituting its radius of localization to Eq. (Q). 
Since at H ~ H c \ we should keep nonlinearity, we can- 
not solve Eq. (Q) analytically. However, MP localization 
radius p can be quite accurately calculated variationally 
with the trial MP wave function in the simplest possible 
form H 



exp 



(15) 



In Van Vleck DMS such variational calculation had 
been done in |b[. For H c i we obtain from dll) 



gip.BH C 20B.{T = 0) = JSpxi 



(14) 



The "universal" (i.e. independent of J and xi) temper- 
ature dependence of H C 20B, is reported in Fig.l. Our 
estimations show that elctron cannot autolocalize also 
in orientational DMS Cdi_ x Mn x Te. Equation ( |l4| ) per- 
mits to estimate H'^or f rom above. Taking into ac- 
count that MP in orientational DMS can be formed at 
xi > 17%, we obtain (for S = 5/2 and xi = 20%) for hole 

H c20R ~ 1800 Tesla - 



H c i — 



87.4 Tesla 



Pi 



Po= P 



Ae 



(16) 



gives po = 0.75 



For xi — 4% variational calculation 
and H^yy = 155 Tesla. 

Thus at xi — 4% vortex lattice for hole MP in Van 
Vleck DMS possibly occurs in wide range of magnetic 
fields 155 Tesla < H < 1300 Tesla. 

We have for H c i in the case of orientational DMS 



H c i — 



6000 Tesla 



Pi 



Po 



(17) 
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For xi = 20% and T = 10K p = 14 and H' c \ OR = 30 
Tesla. 

Thus at xi — 20% and low temperatures vortex lat- 
tice for hole MP in orientational DMS possibly occurs 
for magnetic fields 30 Tesla < H < 1800 Tesla. This in- 
terval is by order of magnitude the same as that in Van 
Vleck DMS. The difference is that in orientational DMS 
we may "adjust" (e.g. getting lower) critical fields by 
temperature. 

We have shown that quantizing magnetic field " forces" 
free spherically symmetric MP in DMS to split into the 
vortex lattice very similar to that in type II supercon- 
ductors. Vortex lattice possibly exists at magnetic fields 
much higher than those for type II superconductors. It 
occurs for Ja < for Van Vleck DMS and for Ja > 
for orientational ones. The values of H c i and H C 2 for 
orientational DMS can be "adjusted" by temperature, 
which is impossible for the Van Vleck DMS. This per- 
mits to hope, that vortex lattice can be detected more 
easily in DMS with orientational paramagnetism. As to 
Van Vleck DMS, in spite of the fact that H c \ and H c2 are 
high, they are accessible by modern experimental equip- 
ment. We think also, that it is possible to find Van Vleck 
DMS, where considered phenomena can occur at lower 
magnetic fields. 
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